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ABSTRACT 

Double-stranded ribonucleic acid-activated protein 
kinase (PKR) downregulates translation as a de- 
fense mechanism against viral infection. In fish 
species, PKZ, a PKR-like protein kinase containing 
left-handed deoxyribonucleic acid (Z-DNA) binding 
domains, performs a similar role in the antiviral re- 
sponse. To understand the role of PKZ in Z-DNA 
recognition and innate immune response, we per- 
formed structural and functional studies of the Z- 
DNA binding domain (Za) of PKZ from Carassius 
auratus (caZa PK z). The 1 .7-A resolution crystal struc- 
ture of caZapi<z:Z-DNA revealed that caZ« PK z shares 
the overall fold with other Za, but has discrete struc- 
tural features that differentiate its DNA binding mode 
from others. Functional analyses of caZ« PK z and its 
mutants revealed that caZa PK z mediates the fastest 
B-to-Z transition of DNA among Za, and the mini- 
mal interaction for Z-DNA recognition is mediated 
by three backbone phosphates and six residues of 
caZapKz- Structure-based mutagenesis and B-to-Z 
transition assays confirmed that Lys56 located in 
the p-wing contributes to its fast B-to-Z transition 
kinetics. Investigation of the DNA binding kinetics 
of caZapKz further revealed that the B-to-Z transition 
rate is positively correlated with the association rate 
constant. Taking these results together, we conclude 
that the positive charge in the (i-wing largely affects 



fast B-to-Z transition activity by enhancing the DNA 
binding rate. 

INTRODUCTION 

The regulation of protein synthesis is one of the most impor- 
tant events under stress conditions. The initiation of trans- 
lation is tightly regulated as the rate limiting phase via the 
phosphorylation of eukaryotic initiation factor 2a (eIF2a) 
(1) by various protein kinases including hemeregulated in- 
hibitor kinase, ribonucleic acid (RNA)-dependent protein 
kinase (PKR), PKR-like endoplasmic reticulum kinase and 
general control nonrepressed 2 (2,3). These kinases sense 
different cellular stress signals such as heme depletion, viral 
infection, Endoplasmic reticulum (ER) stresses and amino 
acid starvation, and ultimately inhibit protein synthesis and 
proliferation (2). Among these kinases, PKR plays an im- 
portant role in the innate immune response against viral in- 
fections by recognizing double-stranded RNA (dsRNA) in 
the cytosol. PKR has a dsRNA binding domain (dsRBD), 
which consists of two dsRNA binding motifs, and an eIF2a 
kinase domain at the N- and C-termini, respectively (4). The 
binding of dsRNA to dsRBD induces the dimerization and 
autophosphorylation of PKR (5,6), and the activated PKR 
leads to apoptosis, inflammation and inhibited translation 
(4). 

Functional orthologs of PKR have been identified in sev- 
eral fish species such as Gobiocypris rarus, Salmo salar, 
Danio rerio, and Carassius auratus (goldfish) (7-10). Since 
these orthologs contain two left-handed deoxyribonucleic 
acid (Z-DNA) binding domains (ZBDs) instead of dsRBD 
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of PKR, they were named as protein kinases containing 
ZBDs (PKZ). Although the sensor domains of PKZ and 
PKR are different, their kinase domains are well con- 
served (8). Similar to PKR, PKZ has been known to reduce 
the overall translational level by phosphorylating Ser51 of 
eIF2a when it is activated by Z-DNA binding (8). There- 
fore, PKZ is a functional ortholog of PKR but it recognizes 
Z-DNA instead of dsRNA. 

Left-handed Z-DNA is a higher energy conformation of 
the double-stranded DNA helix. Unlike canonical right- 
handed B-DNA, which is common in biological DNAs, 
Z-DNA adopts the zigzag arrangement of the phosphate 
backbone as a consequence of alternating anti and syn 
conformations of glycosidic bonds (11). Thus, alternating 
purine and pyrimidine sequences energetically favor the Z- 
DNA formation (12). Z-DNA can be stabilized by reduc- 
ing the electrostatic repulsion between phosphate groups 
in the backbone by high ionic strength, positively charged 
chemicals or modification on bases (13,14). Non-aqueous 
solvents such as ethanol, methanol and ethylene glycol 
also stabilize Z-DNA (15). Besides those physicochemi- 
cal factors, negative supercoiling and Z-DNA binding pro- 
teins (ZBPs) induce Z-DNA under physiological conditions 
(16,17), and the in vivo presence of Z-DNA has been demon- 
strated by Z-DNA-specific antibodies or ZBPs. So far, four 
ZBPs have been identified, each of which contains one or 
two ZBDs (Za): adenosine deaminase acting on RNA-1 
(ADAR1), DNA-dependent activator of interferon (IFN) 
regulatory factors (DAI, also known as DLM-1, and ZBP- 
1), protein kinase containing ZBD from fish (PKZ) and vi- 
ral protein E3L (9,18-20). 

Interestingly, all known ZBPs, ADAR1, DAI, E3L and 
PKZ, are involved in the innate immune response (17,21). 
Mammalian DAI recognizes cytosolic foreign DNA and 
consequently activates the inflammatory signal pathway 
through IFN and NF-kB signaling (22,23). The overexpres- 
sion of ADAR1 or E3L suppresses the DNA-mediated pro- 
duction of IFN in mouse embryonic fibroblast cells (24). 
Hence, the competition between ADAR1 and DAI for cy- 
tosolic DNA has been suggested as a regulatory mechanism 
of the innate immune response. However, the distribution of 
DAI and PKZ is known to vary in different species; that is, 
fish have only PKZ, while mammals have only DAI, sug- 
gesting that PKZ and DAI may perform similar roles as cy- 
tosolic DNA sensors. 

To date, four crystal structures of ZBDs in complex 
with Z-DNA are available: Za domain of human ADAR1 
(Zoladari) (25), Za domain of mouse DAI (Za D Ai) ( 19), Za 
domain of Yaba virus (yabZa E3L ) (26) and Z(3 domain of 
human DAI (hZ(3oAi) (27). Za domains commonly have a 
winged helix-turn-helix (wHTH) DNA binding motif and 
recognize Z-DNA in a conformation-specific manner. Ac- 
cordingly, most of the interactions involve backbone phos- 
phates (17). Despite the similarity in overall structure, the 
binding modes and activities of Za domains are quite dis- 
tinctive (17). For example, it was proposed that the differen- 
tial Z-DNA binding activities of viral Za domains are cor- 
related with pathogenicity (28). 

caZapKz has limited sequence identity with other Za 
domains, 28% for hZaADARi, 20% for hZa D Ai, 23% for 
hZ(3 D M and 22% for yabZa E3L , respectively. However, the 



residues involved in DNA backbone recognition are rel- 
atively well conserved with few exceptions (Figure 1) (9). 
It has already been proven that caZa PK z binds to double- 
stranded (ds) deoxy-hexanucleotide with three repeats of 
cytosine-guanine deoxy-dinucleotide (dCdG) in negatively 
supercoiled DNA and converts ds-(dCdG)3 from B-DNA 
to Z-DNA in vitro (29,30), but little is known about the de- 
tailed Z-DNA binding mode of caZa PK z- Therefore, struc- 
tural studies on Z-DN A recognition are essential to under- 
stand its binding mode to Z-DNA and its relevance to the 
biological function of caZa PK z- In the present study, we de- 
termined the crystal structure of caZa PK z complexed with 
Z-DNA at 1.7-A resolution and observed that caZa PK z me- 
diates B-to-Z transition with the fastest kinetics among the 
known Za domains. Structural analyses combined with the 
functional assay of caZa PK z mutants provided a basis for its 
Z-DNA binding mode and higher rate of B-to-Z transition. 

MATERIALS AND METHODS 

Protein preparation 

The Za domains of human DAI (hZaDAi), human ADAR1 
(hZaADARi), E3L homolog of the yaba-like disease virus 
(yabZa E3L ) and C. auratus PKZ (caZa PK z) were prepared 
as reported previously (19,25,26,30). The detailed proce- 
dures are described in the Supplementary data. 

Preparation of duplex oligonucleotide 

Deoxyoligonucleotides 5'-dTdCdGdCdGdCdG-3' 
[dT(dCdG) 3 ] and 5'-dCdGdCdGdCdGdCdGdCdGdCdG- 
3' [(dCdG) 6 ] were purchased from Integrated DNA 
Technologies (IDT, CA, USA). All oligonucleotides were 
in 50-mM Tris-HCl pH 8.0, 50-mM NaCl and 1-mM 
ethylenediaminetetraacetic acid. Duplex deoxyoligonu- 
cleotides were prepared by annealing an equimolar mixture 
of complementary sequences and dialyzing against buffer 
A. After purification using a MonoQ column (GE Health- 
care, NJ, USA), the concentration of the duplex DNA was 
calculated spectroscopically. 

Circular dichroism measurement 

The B-to-Z transition was monitored by the circular dichro- 
ism (CD) spectrum at 25°C. For each measurement, 7.5 (xM 
of ds-(dCdG)6 in buffer A was used. Wild-type or mutant 
caZa PK z was incubated with ds-(dCdG)6 at final concen- 
trations of 0 |xM (Ox), 7.5 \iM (lx), 15 u,M (2x), 30 |xM 
(4x), 45 \iM (6x) and 60 uM (8x) at 25°C for 1 h. CD 
spectra were recorded using a Jasco J-810 CD spectrometer 
between 230 nm and 320 nm at 1-nm intervals averaged over 
2 s. The kinetics of the B-to-Z DNA transition of caZa P Kz, 
hZaADARi , hZa D Ai, yabZa E3L and their mutants were mon- 
itored by the time course recording of CD spectra at 255 nm 
for 2000 s in the presence of \5-\iM ds-(dCdG) 6 . 

Structure determination 

In our previous study, we confirmed that caZa PK z induces 
the B-to-Z transition of the ds-dT(dCdG) 3 at 2:1 (pro- 
tein:DNA) molar ratio and it forms a stable complex with 
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Figure 1. Overall structure of caZapKz in complex with ds-dT(dCdG)3 and comparison with other Za domains. (A) Sequence alignment of caZapKz 
with human ADAR1 (1iZc<adari). human DAI (hZctDAi) and yaba poxvirus E3L (yabZaE3L). and Z(3 domain of human DAI (hZ^DAi)- The secondary 
structures of caZapicz are drawn at the top of the sequences. The residues interacting with the Z-DNA backbone are marked with filled triangles. (B) 
Overall structure of the caZapKz :Z-DNA complex. Double-stranded DNA and two Za proteins are generated by C2 crystallographic symmetry. Proteins 
are colored slated blue, and backbones and bases of DNA are colored red and gray, respectively. The manganese ions are shown as purple spheres. (C) 
Structural comparisons of the Ca traces of caZapKz and other Za domains. Protein structures of caZapicz, hZaADARl (PDB ID: 1QBJ), yabZc<E3L (PDB 
ID: 1SFU), mZarjAi (PDB ID: 1J75) and IiZPdai (PDB ID: 3EYI) are drawn in red, orange, green, yellow and magenta, respectively. 



Z-DNA at the same molar ratio (30). We crystallized the 
caZotpitz:Z-DNA complex in the C2 space group using 30% 
PEG 1 500 and 1 5-mM MnCl 2 and collected diffraction data 
at 1.7-A resolution. The structure of caZa PK z in complex 
with Z-DNA was determined by the molecular replacement 
method. The initial solution was obtained by MOLREP 
(32) using hZaADARl complexed with ds-dT(dCdG)3 (PDB 
ID: 1QBJ) as a search model. Rigid body refinement of the 
initial solution and density modification were performed 
using REFMAC (33) and DM (34), respectively. Iterative 
model building was done using COOT (35). The final re- 
finement statistics are given in Table 1. The quality of the 
structure was checked using MolProbity (36). All structural 
figures were generated using PYMOL (37). 

RESULTS 

Overall structure of caZa PK z in complex with Z-DNA 

To understand the difference in the B-to-Z transition rates 
of Za domains in relation to their Z-DNA binding modes, 
we determined the crystal structure of caZa PK z in com- 
plex with dsdT(dCdG)3 oligonucleotide at 1.7-A resolu- 



tion with R factor of 17.2% and R[ ree of 22.6%. One 
caZapKz molecule bound to a single-stranded DNA was 
modeled in an asymmetric unit. The current model com- 
prises 62 residues (from Ser2 to Met63) of caZa PK z, seven 
nucleotides of single-stranded dT(dCdG)3 [ss-dT(dCdG)3] 
and one manganese ion. By the C2 symmetry operation, 
two caZapKz molecules bound to both sides of the Z-DNA 
in a head-to-tail orientation were generated (Figure IB). 

The overall structure of caZa PK z is well conserved with 
three (3 strands and three a helices, like other members 
of the Za family (Figure 1). With a wing between the (32 
and (33 ((3-wing), caZa PK z adopts the canonical wHTH 
structure (Figure IB). Structure comparison revealed that 
caZa P Kz is a homolog of other known ZBDs: the root- 
mean-square deviations (RMSDs) from IiZoiadari (PDB 
ID: 1QBJ) (25), mZa DA i (PDB ID: 1J75) (19), yabZa E 3 L 
(PDB ID: 1SFU) (26) and hZ(3 DA i (PDB ID: 3EYI) (27) 
are 0.75 A, 1.57 A, 1.23 A and 0.94 A for 56 Ca atoms, re- 
spectively (Figure 1C). The RMSD versus residue plot also 
indicates that caZa PK z and hZa A DARi have no noticeable 
change throughout the modeled residues. For mZa DA i or 
yabZa E3L , structural deviations from caZa PK z were mainly 
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Table 1. Data collection and refinement statistics of caZapKz/Z-DNA complex 



Native 



Unit cell parameters 
Space group 
Unit cell 



Data collection 
Beamline 

Resolution range (A) 
Redundancy 
Completeness (%) 

-Emerge (%) b 

I/a (I) 

Refinement 

Resolution range (A) 

Number of reflections working sets 

Number of reflections test sets 

Number of protein atoms 

Number of DNA atoms 

Number of waters 

Number of metal ions (manganese) 

«work/«free (%) C 

RMSD bond lengths/angles (A) 
Average B-factor (A 2 ) 
Ramachandran plot (%) d 



C2 

a = 55.37 A 
b = 49.47 A 
c = 29.58 A 
P = 97.22° 

BL4A, PAL 

50.0-1.70(1.73-1.70)* 

3.4 (3.0) 

97.5 (89.9) a 

6.1 (32.5) a 

43.3(4.1) a 

27.5-1.70(1.94-1.70) 

8131 (771) 

382(33) 

504 

139 

70 

1 

17.2/22.6(19.2/25.6) 
0.022/2.226 

25.3(protein)/21.9(DNA)/32.7(water) 
100.0/0.0/0.0 



a The values in the parentheses of the resolution range, completeness, -R mer ge and I/a(I) correspond to the last shell. 

b ^ merg e(/) = XjiklXj'W'^'Oj - A/*W)]l/EhklA«] is the jth measurement of the intensity of reflection hkl and </(hkl)> is the averaged intensity. 

C R = Xjiklll-fobsl - locale II /Hhkll-fobsN where R[ lee is calculated without an s cutoff for a randomly chosen 5% of reflections, which were not used for structure 

refinement, and -R wor k is calculated for the remaining reflections. 

d Percentage of residues in favored/allowed/outlier regions calculated by MolProbity (36). 



observed near the oil helix (Figure 1C and Supplementary 
Figure SI). In contrast, caZa PK z and hZ(3 DAI showed a no- 
ticeable difference in the N-terminal half of a3. The (3-wing 
structure of caZa PK z is also largely different from that of 
mZa DA i and hZ(3 DAI (Supplementary Figure SIB). Struc- 
tural heterogeneity in the al helix is expected since this he- 
lix is least conserved in sequential and structural aspects. 
The structural deviation in the a3 helix and the (3-wing is 
directly related to diversity of Z-DN A binding mode since 
both regions have a major role in recognizing Z-DNA. For 
example, the conformation movement of a3 helix to 3 (10) 
helix has been reported in the crystal and nuclear magnetic 
resonance structures of hZ(3oAi (27,38). 

In the complex structure, dsDNA molecules generated 
by the C2 symmetry operation adopt a canonical left- 
handed Z-DNA conformation (Supplementary Table SI). 
Like other Za:Z-DNA complexes, all three guanines have 
the syn conformation, and the three cytosines have the anti 
conformation. Guanosines 2 and 6 (G2 and G6) adopt 4'- 
exo sugar pucker instead of 3'-endo sugar pucker, as fre- 
quently observed for guanine sugars in Z-DNA (Supple- 
mentary Table S2). All seven nucleotides of ds-dT(dCdG)3 
including 5' deoxythymidine (TO) were modeled (Supple- 
mentary Figure S2A), whereas the 5'-dT overhang was not 
modeled due to weak electron density in the structures of 
other Za:ds-dT(dCdG) 3 complexes. Interestingly, one man- 
ganese ion is coordinated with the backbone phosphate 
(P0), N7 of guanine 2 (G2) and four water molecules with 
an octahedral coordination geometry (Supplementary Fig- 
ure S2B). 



caZapKz has a unique Z-DNA binding mode 

The interaction between caZa PK z and Z-DNA is mediated 
by five residues in the a3 core and four residues in the (3- 
wing. Among them, Lys34, Asn38 and Arg39 in the a3 core 
form direct or water-mediated hydrogen bonds with back- 
bone phosphates of Z-DNA, as observed in the hZctADARi 
structure (25) (Figures 1 and 2). Tyr42, which is known as 
the most critical residue for Z-DNA binding by forming the 
CH-tt stacking with the guanine 4 (G4), is conserved in 
structural and sequential aspects (Figures 1A and 2B). It is 
interesting that caZa PK z has Ser35 in the a3 core, unlike the 
positively charged Lys or Arg in other Za domains, which 
can form electrostatic interaction with the fifth backbone 
phosphate (P5; Figures 1 A and 2). Due to the shorter length 
and lack of positive charge, Ser35 does not form ionic inter- 
actions (Figure 2A and C), but does have hydrogen bond- 
ing with the P3 phosphate. The [3-wing of caZa PK z is com- 
posed of Lys56, Pro57, Pro58 and Trp60. As in other Za do- 
mains, Trp60 forms a hydrophobic interaction with Tyr42 
(19,25,26). The two proline residues, Pro57 and Pro58, in- 
teract with PI and P2 phosphates via hydrophobic interac- 
tions. It is notable that caZa P Kz has Lys56, whereas other 
mammalian Za domains have a polar residue like Ser or Thr 
at the corresponding position. This substitution appears to 
be associated with the unique interaction mode of caZa PK z, 
that is, Lys56 interacts not only with PI but also with P0 
(Figure 2A and C), thereby stabilizing TO and Mn 2+ ion 
(Supplementary Figure S2B). Accordingly, the comparison 
of charge distribution surfaces of Za domains reveals that 
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Figure 2. The Z-DNA binding modes of caZapicz and hZaADARl ■ (A) The binding interface between caZapKz and Z-DNA. The DNA binding interface 
of protein is depicted as a blue ribbon and ball-and-stick models. The residues involved in Z-DNA binding are labeled. The backbone and bases of DNA 
are drawn as red and gray ball-and-stick models. The manganese atoms are represented by magenta spheres. The water molecules are represented by green 
dots. (B) The Z-DNA binding interface of hZaADARl ■ The DNA binding interface of protein is depicted by a green ribbon and ball-and-stick models. The 
residues involved in Z-DNA binding are labeled. The backbone and bases of DNA are drawn as red and gray ball-and-stick models. The water molecules 
are represented by orange dots. The schematic drawings of the DNA binding interfaces of caZapKz (C) and hZaADARl (D) with Z-DNA. The interaction 
residues of caZapKz and hZaADARl are represented in blue (C) and green boxes (D), respectively. The phosphate backbones are numbered from PO to P5 
in red circles. Hydrogen bonds, hydrophobic interactions and it-it interaction are represented by dashed lines, open circles and closed circles, respectively. 
Green (C) and orange (D) ovals stand for water molecules in the caZapKz and hZaADARl structures, respectively. 
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only caZapKz has positively charged wing structure (Sup- 
plementary Figure S3). Current structural interpretations 
imply that the Z-DNA binding mode of caZapKz is distin- 
guished from those of other Za domains, which is due to 
Ser35 and Lys56. Consequently, caZa PK z recognizes PO 
P4 phosphates of Z-DNA, while other Za domains interact 
with P1-P5 phosphates (Figure 2). 

The minimal recognition motif of caZapxz for recognizing 
Z-DNA 

It was proposed from the crystal structures of hZa ADA Ri 
in complex with non-CG repeat of Z-DNA that most con- 
served interactions are mediated by three backbone phos- 
phates, PI, P2 and P3, and five conserved residues, Asnl73, 
Tyrl77, Prol93, Prol92 and Trpl95 (39). It is clear that 
caZapKz seemingly has a different binding mode to inter- 
act with Z-DNA. To elucidate the role of each residue in 
the DNA binding surface and to identify the minimal DNA 
recognition motif in caZa PK z, we mutated the key inter- 
acting residues and monitored the B-to-Z transition ac- 
tivity of mutants at various molar ratios of protein/DNA 
([P]/[N] ratios) (Supplementary Figure S4). Far-ultraviolet 
CD spectra (190-240 nm) confirmed that the overall folding 
of each caZa PK z mutant was not altered, except for Y42A 
and W60A mutants (Supplementary Figure S5). Since CD 
changes at both 255 nm and 292 nm mostly represent the 
B-to-Z transition, we plotted the ellipticity of DNA at 255 
nm and 292 nm in the presence of mutant proteins with 
various [P]/[N] ratios (0-8) (Supplementary Figure S6) to 
compare the B-to-Z transition activity of each mutant. CD 
spectra showed that four molar excess of caZa PK z can fully 
convert ds-(dCdG)6 to Z-DNA, and thus the stoichiome- 
try between caZa PK z and the ds-(dCdG)6 is estimated to 
be 4:1 (Supplementary Figure S4A). Therefore, to simplify 
the comparison, the ellipticity changes at 255 and 292 nm 
from those of B-DNA at the [P]/[N] ratio of 4 were plotted 
for each mutant (Figure 3). In addition, we introduced the 
names of backbone phosphates, which are involved in the 
binding to each mutant in the same plot. Based on their el- 
lipticity change, each mutant was grouped into three. Group 
1 included mutants S35A, K56A and R39A and exhibited 
the largest ellipticity change, suggesting the minimum activ- 
ity change (Figure 3 and Supplementary Figures S4C, E, G 
and S6). The double or triple mutants belonging to group 2 
showed the intermediate displacement, implying that these 
mutants are functionally less active. Mutants in group 3 
appeared to lose their B-to-Z transition activity, which ac- 
counts for the residues essential for the B-to-Z transition. 
For example, the replacement of either Lys34 or Asn38 to 
alanine in group 3 resulted in dramatic reduction of the B- 
to-Z transition (Figure 3 and Supplementary Figure S4B 
and D). In the caZa P Kz:Z-DNA complex, P2 and P3 form 
a wide interaction network with many residues in caZa PK z 
(Figure 2C). Asn38, Tyr42, Pro58 and Trp60 are the bind- 
ing partners of P2, while P3 is recognized by Ser35, Lys34, 
and Asn38 (Figure 2). Accordingly, mutations of residues 
in this interaction network severely damage protein: DNA 
interaction. For example, K34A and N38A showed sig- 
nificantly decreased B-to-Z transition activity (Figure 3). 
These results suggest P2 and P3 interactions are both crucial 
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Figure 3. B-to-Z transition induced by caZapicz mutants. The 
Y-axis and X-axis display the ellipticity (mdeg) difference dur- 
ing B-to-Z transition at 292 nm and 255 nm, respectively. The 
wild-type, K34A, S35A, N38A, R39A, Y42A, K56A, W60A, 
R39A/K56A, S35A/R39A/K56A, K56A/P57A, S35A/K56A/P57A, 
R39A/K56A/P57A and S35A/R39A/K56A/P57A mutants of caZa PK z 
were tested at a [P]/[N] ratio of 4. The plots are clustered in three groups 
by using k-means clustering algorithm. The names of phosphate groups 
inside of parentheses indicate the possible interaction between the residues 
in the mutants and the backbone phosphates. An asterisk indicates the 
phosphate groups that lose one of their residue interactions by mutation. 
In this case, their interactions with proteins partially remain. 



for protein:Z-DNA interaction. This plot also proved that 
Tyr42 and Trp60 have the most critical role in DNA binding 
since the single mutation at those sites almost completely 
abolished the B-to-Z transition function (Figure 3 and Sup- 
plementary Figure S4F and H). Previously, it was suggested 
that these mutants stabilize Z-DNA conformation by form- 
ing a CH-tt interaction (25,26). However, our data suggest 
that they might be also important for the folding of Za do- 
mains since the alanine mutant of Try42 or Trp60 seems to 
have the altered secondary structures compared with wild- 
type protein (Supplementary Figure S5). 

Based on the functional study of multisite mutants, we 
were able to identify the minimal recognition motif of 
caZa PK z and interacting phosphates. The K56A/P57A and 
S35A/K56A/P57A mutants were expected to lose the in- 
teraction with P0 and PI, and showed about 50% reten- 
tion activity relative to the wild-type activity (Figure 3 and 
Supplementary Figure S7). The R39A/K56A/P57A mu- 
tant, which would not interact with P0, PI and P4, had no 
detectable B-to-Z transition activity (Figure 3 and Supple- 
mentary Figure S7). In addition, the R39A/K56A mutant, 
in which the interaction with P0 and P4 is abolished, showed 
about 80% activity (Figure 3 and Supplementary Figure 
S7). These results indicate that P2 and P3 backbone inter- 
actions are not sufficient to stabilize the Z-DNA conformer 
and suggest that interactions with at least three backbone 
phosphates, PI, P2 and P3 or P2, P3 and P4, are required. 
Therefore, the structural geometry of residues involved in 
either PI, P2 and P3 or P2, P3 and P4 recognition is a likely 
minimal structural motif of Z-DNA recognition (Figure 4). 
On the protein side, Lys34, Asn38, Tyr42, Pro58 and Trp60, 
which are involved in P2 and P3 recognition, Pro57 for PI 
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Figure 4. The Z-DNA binding surface of caZapKz- (A) Electrostatic 
distribution of caZapKz- Positively charged, negatively charged and hy- 
drophobic regions are depicted in blue, red and white, respectively. P2 and 
P3 are labeled in red, and PI and P4 are in black. (B) The Z-DNA bind- 
ing surface of caZapKz- Key residues involved in P2 and P3 binding are 
colored orange, and PI or P4 binding residues are shown in yellow. 



binding and Arg39 for P4 binding compose the minimal 
binding motif for Z-DNA binding (Figure 4). 

B-to-Z transition induced by caZa PK z is faster than that in- 
duced by other Za domains 

We examined the B-to-Z transition kinetics by measuring 
the CD spectra at 255 nm for 2000 s. The B-to-Z transi- 
tion rate and half transition time of ds-(dCdG)6 induced 
by caZapKz were calculated to be 15.3 (ms _1 ) and 45.2 s, 
respectively, by fitting the time course CD using non-linear 
regression analyses (Figure 5 and Supplementary Table S3). 
However, B-to-Z transition rates of hZa AD ARi, hZa D Ai and 
yabZa E3L were 5.89, 2.51 and 2.39 (ms -1 ), respectively (Fig- 
ure 5). Accordingly, the half transition times of hZa ADA Ri, 
hZctDAi and yabZa E 3L were 117.6, 275.9 and 290.5 s, re- 
spectively (Supplementary Table S3). Therefore, the B-to- 
Z transition rate of caZa PK z is 2.6 times faster than that 
of hZaADARi, and 6.1 and 6.4 times faster than those of 
hZoiDAi and yabZa E3L , respectively. 

B-to-Z transition activity of caZa PK z mutants 

The time course B-to-Z transition by caZa P Kz mutants 
was monitored by CD spectroscopy in order to provide the 
molecular basis for the high transition rate of caZa PK z- The 
replacement of Ser35 or Arg39 in the a3 core had only a 
marginal effect on the B-to-Z transition rate (Supplemen- 
tary Figure S7 and Supplementary Table S4). However, mu- 
tations in residues in the (3-wing significantly affected the 
B-to-Z transition rate; the calculated half transition time 
was 239 s for the K56A mutant. The activity of the double- 
mutant K56A/P57A was more drastically reduced, result- 
ing in a half transition time exceeding 527 s (Supplemen- 
tary Figure S7 and Supplementary Table S4). These results 
confirmed that the B-to-Z transition rate is largely affected 
by the interaction between the (3-wing of caZa PK z and Z- 
DNA. While residues in (3-wing affect the B-to-Z transition 
rate, K34, N38, Y42 and W60 seem to be essential for B- 
to-Z transition, since their substitution to Ala almost com- 
pletely demolished the transition (Figure 3 and Supplemen- 
tary Figures S4 and S7). We further investigated the roles 
of Ser35 and Lys56 by measuring the transition activity of 
S35K and K56T mutants since their corresponding residues 



A 




Figure 5. B-to-Z transition kinetics of Za domains. (A) CD spectra of 15 
|jlM of ds-(dCdG)6 were monitored at 255 nm in the presence of caZapKz 
(brown), hZaADARi (orange), hZaoAi (yellow) and yabZa£3L (green). (B) 
The calculated B-to-Z transition rates (s _1 ) with error bars from caZapKz 
(brown), hZaADARi (orange), hZa D Ai (yellow) and yabZa£3L (green). 



in hZa ADA Ri are Lys and Thr, respectively, in the sequen- 
tial and structure alignments (Figures 1A and 2). On the 
other hand, we also made hZaADARi mutants by replacing 
Lys 170 and Thrl91 with Ser and Lys, respectively, to make 
caZa P Kz -mimicking mutants. The time course CD spec- 
tra of caZa P Kz mutants (K56T and S35K/K56T), which 
have no positive charge at the [3-wing, showed a decreased 
B-to-Z transition rate (Figure 6A and B, and Supplemen- 
tary Table S5). By contrast, the K56R mutant showed a B- 
to-Z transition rate that was two times faster than that of 
wild-type caZa PK z (Figure 6 and Supplementary Table S3). 
The caZa P Kz-mimicking mutants of hZaADARi, T191K 
and K170S/T191K displayed higher B-to-Z transition ac- 
tivity than the wild-type hZaADARi, while K170S showed 
the reduced activity (Figure 6C and D, and Supplementary 
Table S5), although the difference was not as high as that in 
caZa PK z- These results confirm that charge-charge interac- 
tion between the (3-wing and DNA backbone is important 
for the fast kinetics of Za. Furthermore, we observed that 
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Figure 6. B-to-Z transition activity of caZapKz and hZaADARi ■ (A) Time course CD spectra of caZapKz and S35K, K56T, K56R and S35K/K56T 
mutants. (B) The calculated B-to-Z transition rates (s" 1 ) of the wild-type, S35K, K56T, S35K/K56T and K56R caZapKz are represented by brown, 
orange, lime, green and turquoise colors, respectively. (C) Time course CD spectra of hZc<ADARl and K170S, T191K and K170S/T191K mutants. (D) The 
calculated B-to-Z transition rates (s _1 ) of the wild-type, K170S, T191K and K170S/T191K mutant hZaADARl are represented by brown, orange, lime and 
green colors, respectively. 



the B-to-Z transition activity was decreased when the re- 
action was performed in the presence of NaCl (Figure 7), 
which further demonstrates that the charge-charge interac- 
tion is critical for the B-to-Z transition. The reduction of 
B-to-Z transition rate is more severe in caZapKz than in 
hZctADARi, suggesting that the effect of charge-charge in- 
teraction on B-to-Z transition activity plays a more critical 
role in the case of caZot PK z (Figure 7). 

DNA binding kinetics of ZBPs 

By functional studies of the mutants, we confirmed that the 
positive charge in the (3-wing can enhance the B-to-Z tran- 
sition rate. Given this observation, how does the positive 
charge in the (3-wing affect the B-to-Z transition rate? It 
was proposed that Za bind to B-DNA prior to the tran- 
sition to Z-DNA (40,41), so the association rate of Za to 
DNA possibly contributes to the B-to-Z transition rate. Un- 



der this hypothesis, it is expected that a positively charged 
residue can enhance the transition rate by increasing the as- 
sociation of negatively charged DNA molecules to the pos- 
itively charged protein molecules (42). To further analyze 
the correlation between the association rate and the B-to-Z 
transition rate, the kinetic parameters for the binding of the 
caZapKz, hZa A DARi and their mutants to DNA were cal- 
culated from the Bio-Layer Interferometry measurement of 
the interaction between DNA and Za and their correlation 
with the B-to-Z transition rate (s^ 1 ) was examined (Supple- 
mentary Figures S8, S9 and S10; Supplementary Table S6). 
We found that B-to-Z transition rates of the wild-type and 
mutants caZa PK z are proportional to the association rate 
constant k on (Supplementary Figure S9E) that is also pro- 
portional to the number of positively charged residues in the 
Z-DNA binding surface. A similar correlation is also found 
in hZotADARi, although the correlation coefficient (r) is rel- 
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■ B-DNA 
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Figure 7. Salt effect on the B-to-Z transition. (A) Time course CD 
spectra of Z-DNA:caZapKz in various NaCl concentrations (10, 40, 80, 
100, 150, 200 and 250 mM). The ellipticities (mdeg) monitored at 255 
nm are indicated by different colors. (B) Time course CD spectra of Z- 
DNA:hZoiADARi i n the same conditions with the same color schemes as in 
Figure 7A. 



atively lower than that of caZa PK z (Supplementary Figure 
S10 for hZctADARi)- 



DISCUSSION 

We have compared the Z-DNA binding structure and B-to- 
Z transition activity of caZapKz with those of other Za do- 
mains. A structural comparison revealed not only the con- 
served overall binding mode of caZa PK z but also some dis- 
tinctive features that appear to be important to its higher B- 
to-Z transition rate. Unlike other ZBDs, caZa PK z harbors a 
positively charged residue in the (3-wing, Lys56, which plays 
a key role in DNA binding by forming a wide interaction 



with DNA, including bonds with both PI and P0 (Figure 
2A and C, and Supplementary Figure S2B). 

Although the overall charge distribution is not conserved 
among ZBDs, positively charged residues and hydropho- 
bic residues are commonly observed in the ct3 core and 
in the (3-wing, respectively, which play essential roles in Z- 
DNA binding. However, it seems that the positively charged 
residues have a more significant effect on DNA binding 
and the B-to-Z transition than hydrophobic residues. There 
are three basic residues involved in Z-DNA recognition in 
caZapKz (Lys34, Arg38 and Lys56) and hZa A DARi (Lysl69, 
Lysl70 and Argl74) (Figure 1A). However, hZa D Ai and 
yabZa E 3L have only two positive residues in the Z-DNA 
binding surface (Figure 1A). Accordingly, hZa A DARi in- 
duces 2.3 and 2.5 times faster B-to-Z transition than 
hZotDAi and yabZa E3L , respectively (Figure 5B). In con- 
trast, the B-to-Z transition rate of caZa PK z, which has 
one more basic residue in the (3-wing than hZa D Ai and 
yabZa£3L, is 6.1 and 6.4 times faster than that of hZctDAi 
and yabZa E3L , respectively (Figure 5B and Supplementary 
Table S3). These results suggest that the contribution of 
positively charged residues during the B-to-Z transition is 
more significant when they are present in the [3-wing than 
in the a3 helix. Consistently, the K56T mutation in caZa P Kz 
significantly lowered the transition rate from 14.5 ms~' to 
4.68 ms" 1 (Figure 6A and B, and Supplementary Table S5), 
while the addition of a positively charged residue in the (3- 
wing enhanced the activity of hZaADARi from 5.32 ms -1 
to 8.83 ms" 1 (Figure 6C and D, and Supplementary Table 
S5). However, the effect of charged residues in the a3 he- 
lix seems to be not as significant as the effect of charged 
residues in the [3-wing, although the positive charge in the 
a3 helix still contributes to an increased B-to-Z transition 
rate. Furthermore, the S35K mutant showed a higher transi- 
tion rate than the wild-type caZa PK z (Figure 6A and B, and 
Supplementary Table S5), and the K170S mutant showed 
a slightly decreased B-to-Z transition rate compared to the 
wild-type hZa A DARi (Figure 6C and D, and Supplementary 
Table S5). These results are also consistent with the kinetic 
study of the viral Za protein; the positively charged amino 
acids in the (3-wing of Za E3L enhanced the B-to-Z transition 
activity (28). Taking these results together, we conclude that 
the positively charged residue in the [3-wing could facilitate 
the B-to-Z transition (Figure 6). The fastest B-to-Z transi- 
tion rate of the K56R mutant of caZa PK z among wild-type 
and all mutant ZBDs (Figures 5 and 6, and Supplementary 
Table S3) and the reduction of the transition rate under high 
salt conditions (Figure 7) also support our conclusion. 

From a structural point of view, Arg56 in the model of 
the K56R mutant of caZa P Kz can fit nicely into the pocket 
between P0 and PI, suggesting that Arg can form a tight 
interaction with P0 and PI groups (Supplementary Figure 
Sll). From the further analyses of their DNA binding ki- 
netics (Supplementary Figures S8-S 1 0), we propose that the 
increased transition activity is achieved by the enhanced as- 
sociation of ZBDs to DN A through their charge-charge in- 
teraction. However, this correlation was not found when dif- 
ferent Za domains were compared (Supplementary Figure 
S8), which is possibly explained by the fact that the con- 
tribution of a positive charge at a3 and the (3-wing during 
the B-to-Z transition is not the same as that shown in Fig- 
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ure 6 and in Supplementary Table S5. Given these results, 
we conclude that the association rate constant is affected 
by the presence of a positively charged residue on the DNA 
binding surface, which eventually enhances the B-to-Z tran- 
sition rate. Moreover, the contribution of a positive charge 
in the f$-wing is larger than that in the helix a3. 

In terms of structure, unlike the other crystal structures 
of ZBDs complexed with Z-DNA (19,25,31), Mn 2+ ion is 
bound to the caZa PK z:Z-DNA complex. It was proposed 
that nickel and cobalt ions can bind to guanine N7 due to 
the increased solvent accessibility in the Z-DNA conforma- 
tion (43,44) and the binding of manganese ion to Z-DNA 
was recently reported (45). However, among all N7 atoms of 
guanine in caZa PK z, N7 of G2 appeared to specifically par- 
ticipate in the coordination of the manganese ion. The man- 
ganese ion is also coordinated by O IP of CI, which is sta- 
bilized by Lys56, so Za binding is thought to be necessary 
for coordinating the manganese ion. To investigate the role 
of manganese ion in the biochemical function of caZapKz, 
the B-to-Z transition rate of caZa PK z was measured in the 
presence of 15-mM MnCT. However, the manganese ion 
seems to not contribute to the fast B-to-Z transition ki- 
netics of caZotpKz> since the transition rate decreased from 
15.3 (ms -1 ) to 9.49 (ms _1 ). hZotADARi and hZa D Ai showed 
marginally increased B-to-Z transition kinetics in the pres- 
ence of manganese ion (Supplementary Figure S12). There- 
fore, MnCh does not seem to be correlated with the fast B- 
to-Z transition kinetics of caZa PK z, and thus the binding 
of manganese ion to caZa PK z may be considered as a crys- 
tallization artifact since the crystallization buffer contained 
15-mM MnCl 2 . 

It has been postulated that the interaction of ZBP to Z- 
DNA is critical for various biological functions, and thus 
the B-to-Z transition activity of the ZBPs is relevant to the 
biological functions. DAI binding to the cytosolic dsDNA 
turns on the innate immune signal by triggering the acti- 
vation of NF-kB signaling and expression of type-1 IFN 
(23). However, PKZ works as an innate immune sensor by 
recognizing poly-d(GC) as Z-DNA and turning off gen- 
eral translation machinery by phosphorylating eukaryotic 
translation initiation factor 2a (eIF2a) (8,29). Interestingly, 
in fish species, both PKZ and PKR are known to play an 
antiviral role as IFN-inducible eIF2a kinases (46). There- 
fore, we are intrigued by the functional relationship between 
PKZ and PKR in terms of the antiviral response. A re- 
cent study of PKZ and PKR in C. auratus revealed that 
both proteins work independently, but play a cooperative 
role in IFN-mediated antiviral response, although antiviral 
ability offish PKZ was weaker than that of fish PKR (21). 
However, the pull-down experiment with poly I:C, a dsRNA 
mimic, showed that the DNA binding domain of PKZ does 
not bind to poly I:C, while PKR does bind (21). These re- 
sults suggest that PKR and PKZ might independently con- 
tribute to the innate immune response by recognizing differ- 
ent nucleotide substrates in the initial stage of the innate im- 
mune response. Indeed, PKZ-dependent phosphorylation 
of eIF2a and the subsequent translational shutdown are 
largely induced not by dsRNA but by Z-DNA (8,29). Since 
the Z-DNA-induced innate immune responses conducted 
by DAI could be regulated by competition with other ZBPs 
such as ADAR1 or E3L, differential Z-DNA binding activ- 



ity or the B-to-Z transition rate of Za domains in the innate 
immune sensors might be a key to controlling the innate im- 
mune response. Similarly, since fish virus also harbor the 
ZBPs (47), it is also expected that PKZ must compete with 
viral ZBPs to play a role in fish innate immune response 
during viral infection. Therefore, the fast Z-DNA convert- 
ing activity of PKZ might be necessary for competing with 
viral ZBPs. Further studies on the B-to-Z transition activity 
of Za domains in cellular environments will clearly define 
the biological relevance of the structural and kinetic differ- 
ences between caZa PK z and other Za domains. In addition, 
structural and functional studies on intact ZBPs are also re- 
quired for better understanding of the Z-DNA-dependent 
innate immune response. 

The current study also provides an insight into a novel 
ZBD. Z-DNA is known to have various roles in cellular 
events. Z-DNA formation is necessary for opening chro- 
matin structure and activating transcription of target genes 
(48,49). Z-DNA is also known to cause genomic instabil- 
ity (50). For the regulation of diverse functions of Z-DNA, 
more ZBPs with novel ZBDs are expected to mediate vari- 
ous cellular processes associated with Z-DNA. The results 
from this study combined with those from structural and 
functional studies of Za A DARi confirm that Z-DNA can be 
recognized by Za proteins through three backbone phos- 
phates, but the interaction modes for caZa PK z and Za A DARi 
are not the same in detail (39). Considering the size and 
charge distribution of the minimum unit of Z-DNA for pro- 
tein binding and the diversity of the Z-DNA binding mode, 
it is expected that more ZBPs other than the Za domain 
may be present in cells. Consistently, this possibility is par- 
tially proven by identifying putative ZBPs that interact with 
the Z-form of ADAM-12 NRE sequence (51). Therefore, 
the current study not only reveals a structural requirement 
for fast B-to-Z transition but also provides insight into the 
presence of novel ZBPs. 
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